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Abstract 

Hadronization corrections to the thrust and related event shape distributions in the two-jet 
kinematical region of e + e~ annihilation are summarized by nonperturbative shape functions. The 
moments of shape functions are given by universal matrix elements in QCD, which describe the 
energy flow in QCD final states. We show how the nonperturbative structure of these matrix ele- 
ments may be inferred from resummed perturbation theory and Lorentz invariance. This analysis 
suggests the same functional forms for the shape functions as were found in phenomenological 
studies, and sheds light on the physical significance of the parameters that characterize these 
functions. 
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1. Introduction. The predictive power of perturbative QCD is based on the factorization 
of long- and short-distance dynamics, made possible by their quantum-mechanical incoherence. 
The hard scattering of quarks and gluons can be computed in ordinary perturbation theory, while 
soft-scale physics is encoded into nonperturbative functions, which are fit to data, or evaluated 
using models. In a rather limited number of inclusive processes, which admit the operator product 
expansion (OPE), these functions are described by universal inclusive distributions, for example, 
leading-twist parton densities in deeply inelastic scattering. The latter can be defined in QCD in 
terms of correlators of quark and gluon fields on the light cone, and can, in principle, be computed 
using nonperturbative methods. On the other hand, there is a large variety of processes that are 
not completely inclusive, but are nevertheless infrared safe |TJ], and thus calculable in perturbation 
theory. Important examples are event shape distributions in e + e~ final states f2[. Nonperturbative 
physics enters these cross sections as corrections, suppressed by powers of the large energy scale, 
Q. 

Since event shape distributions are weighted, rather than fully inclusive, cross sections, their 
power corrections are not constrained by the OPE, and are generally less suppressed. For instance, 
the mean value of the thrust, T, acquires a nonperturbative correction (£)np ~ A/Q with t = 1— T, 
and A a new QCD scale ||. We may think of this correction as the first term in an expansion 
in 1/Q. The situation is even more complex for differential event shape distributions, such as 
da/dt, in the two-jet kinematic region, t — *■ 0. In this case, one has to deal with an infinite 
number of scales, (da/dt) NP ~ J2k Ak/ (tQ) k . The organization of such a series is equivalent to the 
introduction of a new nonperturbative distribution. These distributions are the shape functions 
||, [|. They provide a successful phenomenological description of differential event shapes, down 
to the two-jet limit, and over a wide interval of energies, as shown in Refs. ^ ffl]. In this paper, 
we show how the functional form found in these phenomenological studies may be motivated from 
QCD. 

A perturbative analysis of the event shape distributions in the two-jet region shows that the 
leading power corrections, l/(tQ) h , are associated with multiple soft, wide-angle gluon emission 
at the momentum scale tQ [Q, Efl. The influence of collinear splittings of quarks and gluons is 
suppressed by an extra power of Q. The leading behavior of the cross section is thus independent 
of the dynamics that produces the internal structure of jets, and is well described by a physical 
picture of the perturbative final state in which two fast moving back-to-back quarks propagate 
through a cloud of soft gluons, behaving as classical sources of color charge. This is equivalent to 
the eikonal approximation for the quark-antiquark pair. 

In eikonal approximation, multiple soft-gluon emission from quarks is described by nonabelian 
phase operators 4> nR [oo,y] and ^ [oo,y]. Here @ n [oo,y] = P exp (ig J °° dan' 1 A^an + y)) is a 
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Wilson line, with its light-like direction defined by a quark momentum, pi = n1JQ 2 /2 or 



P2 = n RyQ 2 /'2, with Q = (P1+P2) the total center-of-mass energy. In a frame where the two-jet 
axis points in the ^-direction, we choose the four- velocities n R = ^=(1, 0, 1) and n^ L = -^(1, 0, — 1), 
with n\ = n 2 R = and tll ■ ur = 1. Combining the eikonal phases of quark and antiquark, we 
obtain 

U(y; A) = T [<P nR [00, y\& nL [00, y] } = T {<P nR [00, y\Q nL [y, 00]} , (1) 

where we have used the unitarity of the phase operators. The symbol T stands for time ordering 
of the gauge fields in the product. The explicit expression for U can then be written as 



U(y; A) = TPexp [ig [ dz^A^y + 



(2) 



where P orders the color indices of gluon fields along the integration path C e + e -: z^ir) = 
t n R 9{—t) + rn^ 6{r) with —00 < r < 00. 

A characteristic of weighted cross sections like the thrust and the heavy jet mass is that they 
distinguish between gluons moving into the right and left hemispheres, defined with respect to 
the plane orthogonal to the two-jet axis (in the overall center-of-mass frame). When the quarks 
are treated in the eikonal approximation, the shape function that describes energy flow for these 
final states is M, |5H 



f(e R , e L ) = E mU(0)\h, • • • , k N )\ 2 S (e R - £ k i>+ ) 5 f 

N \ i£R / \ 



£ l-J2 h- ■ (3) 

ieL J 

Here, the summation goes over soft particles in the final states, with momenta k\, . . . , k^. The 
variables in the shape function, Sr and El, are projections of the total momenta of soft particles 
flowing into the right and left hemispheres onto the corresponding jet directions, 

k + = k ■ tir, k_ = k ■ ul . (4) 

The operator U(0) in Eq. (H) represents radiation from the quark- antiquark pair, in the eikonal 
approximation discussed above. 

To derive a practical operator interpretation for the event shape function (^) in QCD, it is 
convenient to eliminate the summation over final states. To this end, we introduce an energy 
flow operator, £(n), which acts on asymptotic states by measuring the differential distribution of 
particles with energy k° moving in the solid angle dn5(n 2 — 1) = dcosxdf || |9|, 0. We define 



this operator by its action on a final state consisting of iV particles with momenta ki, . . . , k^, 
according to 

JV 

8{n)\k u ...,k N ) = 5^A;°<J(cos0j - cosx)5(<fr j - <p)\h, ...,k N ), (5) 



Figure 1: The pictorial definition of the energy flow operator S(n). The unit vectors n R and 
Hl indicate the directions of the two outgoing jets. The shaded plane orthogonal to these unit 
vectors goes through the annihilation point, and separates the left and right hemispheres. The 
energy flux through the infinitesimal surface element du = dn5 (n 2 — 1) is given by 8{n)duj. 

—* — * 

where spherical angles \9j,<j>j) define the position of the unit vector Hj = kj/\kj\ with respect to 
the jet axis n R , and (x, <fi) define those of n. Introducing a light-like vector n M = (l,n), we can 
represent the total momentum flow through the right (n 3 > 0) and left (n 3 < 0) hemispheres of 
infinite radius in terms of 

V» L = J dnn»6{±n 3 )5 (n 2 - l) E{n) , (6) 

so that 

Pr|/ci, ...,k N ) = Y^ kf\h, ...,k N ), (7) 

j&R 

and similarly for V^. Notice that the total momentum of the final state is given by the sum 
= J dx Ofl (t, x) = Vji + Vl, with the energy-momentum tensor. 
Using these definitions, the shape function, Eq. (^), can be reexpressed as 

f(e R , e L ) = (0\tf(0)5 (e R - £ R ) 6 (e L - £ L ) C/(0)|0) , (8) 

where 

£r,l = J dn5 (n 2 - w RjL (n)£(n) , w RiL (n) = (1 - | cosx|) 9 (±cos%) . (9) 

Defined in this way, the shape function becomes a new QCD distribution, which governs nonper- 
turbative power corrections to a number of differential cross sections in the two-jet region. 

A typical application of the shape function, which illustrates the factorization of soft gluon 
emission from perturbative dynamics, is to the distribution for the heavy jet mass. The heavy jet 
mass is defined by g = (l/<5 2 )max(m|,, m 2 L ), with m R (m^) the invariant mass of the particles 
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moving into the right (left) hemisphere. It may be written in the convolution form 



Qfg (QQ, QQ) ^j T (0) + I def e (e, gQ) — -— ^ ^ , (10) 



cr tot Jo d£ 

where daj T /dg is the resummed perturbative single-jet cross section pTT[| , 7£j T is the corresponding 
'radiator function', defined by dTZ^ T (g)/dg = daj T /dg, and 



/,(e, eg) = 2 jT° cfe'/ (e, e') K p / (q-^\- ( n ) 



The thrust and the C-parameter distributions are given by similar expressions 0, where f g is 
replaced by a single-variable function, related to /(£r,£x) byf] 

/■oo 

/ t (e) = / tfeatfei, /(ea, e L ) 6(e-e L - e R ) = (5 (e - S R - £ L )) . (12) 
Jo 

Corrections to Eq. (|I0"D are suppressed by powers of Q, as discussed above. 

In Ref. M, it was found that a successful description of the differential distributions for thrust, 
heavy jet mass, and C-parameter over a wide range of energies is given by the single shape function, 

t( v JV(a,b) (SRELy- 1 ( e R + e 2 L + 2be R e L \ 

f(e R ,e L ) = (^—J exp j , (13) 

when its parameters are chosen as: a = 2, b = —0.4 and A = 0.55 GeV. The normalization 
factor J\f(a,b) is chosen so that / ds R ds R , El) = 1. Recently, a thorough investigation of 
the differential distribution for thrust |7j] reached similar conclusions in terms of the function f t , 
Eq. (0). In the following, we shall argue that the functional form, Eq. (|13"D, follows from a 
perturbative analysis of moments of the matrix element (|8|), supplemented by Lorentz invariance 
and standard treatment of integrals of the strong coupling over low momentum scales. 

2. Energy correlators. In a manner analogous to the relation between the light-cone ex- 
pansion and moments of deep-inelastic scattering structure functions, the moments of event shape 
functions give matrix elements of certain (in this case, nonlocal) operators in QCD. Specifically, 
from the definition, Eq. (g), the moments of f (er, £_l) produce weighted Green functions of prod- 
ucts of the operators £, inserted into correlators of Wilson lines, 

POO 

/ ds R de L e n r e^ f(s R , e l ) = (£»£?) (14) 
Jo 

. N+M N M 

= y II d ™i 5 (n 2 j - x ) II w r{^ti) II w L {n m )Q{n u n N+M ) , 



where 



<?(ni, . . . ,n N ) = (S(tii) ...£(n N )) (15) 
= ^(0|Tr {t/t (0; A<->) £{n x ) . . . £(n N )U (0; #') } |0> 



2 In what follows we use a convention (Q\W . . .U\0) = 



4 



measures correlations between energy flows in the directions of specified unit vectors fij. Rather 
than study the shape function directly at first, we begin with these, even more general matrix 
elements, from which the shape function (H) may be derived ||. 

We shall analyze below the energy correlators, G{n\, . . . ,n^), from a general point of view 
that has been widely applied to infrared safe quantities, using a variety of related prescriptions 
||. This may be summarized very briefly as follows. Starting with low-order diagrams for the 
quantity in question, we absorb higher-order corrections into the scale of the running coupling. In 
the resulting expression, the integral over the argument of the running coupling from zero up to 
a factorization scale, A, is replaced by a phenomenological parameter. Other, overall kinematic 
integrals, upon which the running coupling does not depend, determine properties like angular 
and rapidity dependence. We will show that this procedure, applied at order a s and a 2 to the 
energy correlators, is adequate to derive the functional form Eq. ( pT3|) for the corresponding shape 
functions, and to give physical interpretations to its dimensionless parameters a and b. 

The perturbative expansions of the correlators of Eq. ([TJ) are illustrated diagrammatically by 
Fig. |^. As the figure shows, these correlation functions correspond to cut diagrams, rather than 
time-ordered Green functions, since the gauge fields entering the eikonal phases U and IP are time 
and anti-time ordered, respectively. To express this fact, we adopt the Keldysh notation, assigning 
"plus" and "minus" superscripts to the fields entering the amplitude (A^) and its conjugated 
counterpart (A^). The propagator of the gauge fields A^ is defined as 

J d 4 xe^{0\A^(x)A^ (0)|0> = ^(q)D^\q), (16) 

with a, (3 = ±, and with the conventional free gluon polarization tensor in covariant gauge, 
Qnv{q) = 9»v - (1 - QqnQv/q 2 - in tnese terms, we have = -[£)("%)]* = (q 2 + iO)' 1 , 

D^ + \q) = D( + ~\—q) = —2iri9( y q°)S(q 2 ). We choose below the Feynman gauge, £ = 1, for 
calculations (although our results are gauge invariant). 

3. The single-gluon approximation. We begin our study of the properties of the correlation 
functions with the lowest-order perturbative expression for the single energy correlator, Q{n). 
Using the definition (|15D and expanding (0) in powers of the gauge coupling, we have 

G{n) =g 2 [ dz» f dz^ (0| A^( Zl )£(n)Ai+\z 2 ) |0> + O (g A ) . (17) 

Going to the momentum representation, and performing the integration along the contour, we 
find 

C j?cx f dk f dk^ f 
G{n) = 2 ^ 2 S J -j— J -p- J d(p k k 5(cos6 k - cos x)S{4> k - f) ■ (18) 

In this formula, we recognize the lowest-order bremsstrahlung spectrum of soft gluons, weighted 
by an additional energy flow factor, according to Eq. (|5]). Using the definitions of the light-cone 
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e^e" 



Figure 2: Diagrammatic representation of the correlation function (|T|). Here the dashed line rep- 
resents the unitary cut and the final states are weighted with the product of factors £{n\) . . . £{n^) 
defined in Eq. (||). 



variables k± = ^= (k ± k%), along with k = v^sin 1 9k, k 3 
d cos 9k/ sin 2 9k- The integration over k + then yields the result 

1 1 



k 2 cot 9k, we find dk + /kj l 



G{n) 



27r sin x 



dk v fe Ppt ( 



(19) 



where we have introduced a resummed perturbative density, ppt (A; 2 ) = ^-f "^2 ^ 0, which mea- 
sures the number of particles produced per unit interval in rapidity and transverse momentum. 
The expression (|T^) exhibits a third-order pole as x ~> or x n - 111 both gluon propa- 

gates in the final state close to the direction of an outgoing jet, and the resulting singularities are 
of collinear origin. 

The simplicity of the x-dependence in Eq. (|TJ^) reflects the underlying boost-invariance of the 
nonabelian phase operators U. All dynamical information is contained in a boost-invariant particle 
density, independent of the rapidity, 



2 k. 



In tan 



X 



(20) 



The resummation incorporated into the running coupling in the step from Eq. (|18f) to Eq. (|19f) is 
the only choice consistent with the boost properties of the correlators. This resummation, which 
summarizes higher-order scale-fixing corrections, a s (j3oa s ) N , is normally justified diagrammati- 
cally, and is used as a starting point for infrared renormalon-inspired arguments . We shall not 
employ the form of the perturbative running coupling, however, nor will we make a specific ansatz 
for a nonperturbative generalization of the integrals in Eq. ([191) . Following Refs. fU, however, 
we replace the resummed perturbative density by a phenomenological one, pnp(& 2 ), below some 
cutoff, A 2 , in the k 2 integral. Whatever the choice of pnp, h should, like its perturbative counter- 
part, be boost invariant. For any such choice of the nonperturbative density, the single-particle 



6 



energy correlator, Q(n), defines a leading-power correction to the energy-energy correlation in 
e + e~-annihilation |12[ , 

(EEC( X ))np = ^G(n) , (21) 



with the rapidity dependence specified by Eqs. ([19]) and (^3), and any choice also produces a 
uniform particle number distribution, d 2 N / drjdk 2 = p(k 2 ), in terms of the rapidity. 

We have thus, starting with the energy correlators, arrived at a "tube" model, as reviewed in 
Ref. ]13[, in which particles are produced with a constant, nonperturbative density at all rapidities. 



The perturbative calculation therefore suggests a realistic dependence on the rapidity, in terms 
of an integral of the density p over nonperturbative scales. This integral, of course, is to be 
replaced by a set of phenomenological parameters. We now use this approach in the computation 
of multiple energy correlators. Perhaps surprisingly, the lowest-order model offers insights here as 
well. 

To the lowest order in perturbation theory, the multiple energy flow correlator Q{n\, . . . ,nj) 
is given by 

g(n 1 ,...,n j )=g 2 [ dz» [ dz» (0| A^(z 1 )£(n 1 ) . . . £{n 3 )A^\z 2 ) |0> + O (g 4 ) , (22) 
whose evaluation leads to 

Q(ni,...,ftj) = ^-^7^2 — Il 5 ( cos Xi - cosxi)<%i - <pi) f dk 2 (Vk?) 3 p PT (fc 2 ) , (23) 

ZiTX Sill ■ o J 



where (xiifi) are the spherical coordinates of the unit vector f^. Notice that this correlator 
vanishes unless all vectors are aligned. Higher order corrections, discussed below, will relax this 
condition. 

Starting from Eq. (|23|) , we easily calculate the total flow into the right or left hemispheres. 
According to (0), it is given by 

o r y[ 2 »7 n AN 

={£%) = I * 2 (^) "-(fc 2 )=2C^A_, (24) 

where in the last relation we give the result for fixed coupling. The vanishing of inter-hemisphere 
correlations, 

= + O (a 2 ) , (25) 

implies that the shape function, /(e^,^), is factorized into the product of two functions, each 
depending on a single energy variable 

/onc-gluon(^i?, £l) = fi^Rjfi^L) , (26) 
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with the function defined as f(£i) = (5(ei — Si)) for i = R,L. Nonfactorizable corrections to 
Eq. (|26|) begin with a^-contributions to multiple energy flow correlators G(ni, . . . ,flj) associated 
with final states with two real gluons. Phenomenological analyses || [7J have shown that power 
corrections to certain event shapes, in particular the heavy-jet mass, are very sensitive to these 
corrections. 

4. Order a 2 and correlations between hemispheres. To study correlations between 
hemispheres at lowest order, a 2 , in G{ni,n 2 ), we restrict ourselves to Feynman diagrams in Fig. 
|2] involving two-gluon cuts. According to the color structure, one can separate abelian and non- 
abelian contributions, G = Ga + Gna- 

The abelian part, GA(ni,n 2 ), arises from QED-like diagrams. Due to the exponentiation 
properties of eikonal lines ||14|| , their contribution is factorized into the product of independent 
single-gluon emissions, 

n 2 C 2 M 2 dh 2 dh 2 

Ga{^u 2 ) = ^ 1 . 3^.3 / ^L^L = G{n l )G{n 2 ), (27) 
47r 4 simxi sin 6 %2 Jo Jtf ~ ' 



where G{n) is given by Eq. (^8[) and where we impose a cutoff, A, on gluon transverse momenta, as 
above. The overall factor (sin \i sin X2)~ 3 corresponds to two independent collinear enhancements, 
when either gluon is emitted nearly parallel to one of the eikonal directions. For now, we study 
Ga^i-i^) at fixed a s . We will comment afterward on the role of the infrared behavior of the 
running coupling. 

The nonabelian part, Gna^i, n 2 ), receives contributions from the crossed ladder, self-energy 
and triple-gluon vertex diagrams (see, e.g., Fig. 1 in Ref. |15|]). The total expression is a somewhat 
complicated function depending on the spherical angles (xi, <fi) and (X2, ¥2)- Closer examination 
reveals, however, that Gna^i,^) depends, up to a prefactor, only on two specific combinations 
of these angles 

oi 2 CfCa ( \ 

GNA(nx,n 2 ) = ~4 . 3 F A 3 — A 2 J r ( cosh(rj 2 - rji),cos(ip 2 ~ fi)) , (28) 
7r z sin Xi sin X2 V / 

with 

1 ~x_\ X2 1 X2 Xi 

cosh(?? 2 — rii) — - tan — cot 1 — tan — cot — , (29) 

w ' 1 2 2 2 2 2 2 y ' 



where the rapidity rji is defined in Eq. (P0f) . This property is a consequence of boost invariance 
along the z-direction, and rotational invariance around it. The explicit expression for the dimen- 
sionless function JF is rather cumbersome, and in fact we do not need it for our purposes. Indeed, 



the correlator G enters into the expression for the moments of the shape function (fT3| ) under the 
integral over azimuthal angles Since the weights w^l (0) do not depend on the (fi, we may 
integrate over the azimuthal angles, to arrive at 

^ ri* a 2 C F C A 

7r 2 sin 3 Xi sin 3 X2 



f d<pi f dip 2 GNA{ni,n 2 ) = , 3 CfC . A 3 — A 2 T{\i) 2 - 771 1 ) , (30) 
Jo Jo n z sin Yi sin y 2 
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where 



{2A) 2 h ^4^4 D^klkl) \ V i' ' ' fcf + fcisinhr/ 

with the function Z) defined as 



h 2 b 2 

1 2 



D ( Jfej, fc|) = 1 + 2 / 2 2 cosh 77 . (32) 



Integrating over transverse momenta, we find 

cosh 77 



J-{rj) = (coth 77 — 1) (2 + 77 cosh 77 coth 7/) H 3 — (277 cosh 2 7/ — 3t? + sinh 

2 sinh 77 v 



- ( cothr? - 1 + — — ) ln(2 + 2 cosh 77) . (33) 
V smh 77 y 

The function JF depends only on the relative rapidity, 77 = \t)i — 772 1 , of the emitted gluons. We 
find that at large 77 it decreases exponentially 

F(rj) = 277e- 3 " + O (e" 4 ") , (34) 

while at small 77 it has a pole, 

^) = ^{§-41n2} + O(77 ). (35) 

These properties may be understood as follows. 

When Xk ~~ > 0, tt, the fc-th gluon is emitted collinear to one of the eikonal directions. At 
first sight, the nonabelian contribution, Eq. (^8) has the same double collinear enhancement as 
the abelian term, fl27|) from the factor (sin%i sinx 2 )~ 3 - The large-77 behavior exp (— 377) in F(r]), 
however, prevents enhancements when the directions of the gluons are widely separated in rapidity. 
The nonabelian term therefore has at most a single collinear enhancement from configurations 
where both gluons approach the direction of one of the eikonal lines. This result is consistent 
with the exponentiation of collinear logarithms in inclusive cross sections, which requires that the 
double-logarithmic collinear enhancement at two loops be given by the square of the enhancement 
at one loop. 

At small rapidity intervals, 1 772 — ??i | — ^ 0, at fixed 77^ 7^ 0,?r, the gluons propagate in the 
same direction. This gives rise to another collinear pole, I/77, from (RST), which is compensated 
by a virtual 0(a s ) correction to single gluon emission. Once combined, they define a distribution 
regularized via a 'plus'-prescription, and give no new collinear enhancement to any infrared safe 
cross section. 
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Let us examine the contribution of the correlator Q(ni,n2) to the total energy flow into the 
left and/or right hemispheres, Eqs. (^4j) and (|25|) . We find that (|24l) receives an 0(a 2 s ) correction, 
while (^5]) develops a non-zero value. The leading order af-contribution to the correlator (£r£l) 
is naturally decomposed into an abelian (~ Cp) and nonabelian (~ CfCa) part. Due to the 
factorization property, (p7[), the former is reduced to the product of two factors, {£r){£l)i eacn 
given by the lowest order expression (|24D . The nonabelian piece defines the irreducible part of the 
correlator (£ R £ L ), 

((£r£l)) = (£r£l)-(£r)(£l)- (36) 

It is given by the real-gluon contribution to Gna^i,^) alone, evaluated in Eq. (|30J). Because 
ni and n 2 belong to different hemispheres, the virtual gluon contribution is non-vanishing only 
for xi — X2 = 7r/2. This corresponds to a single point in the phase space, < xi < 7r /2 and 
< X2 < 7t? and vanishes upon integration over xi,2- Substituting the result of the perturbative 
computation of Eqs. fl3"0|) and (^), one finds after some algebra 

((£ R £ L )) = — 2 C F C A A 2 \ dririe-Wtn) w -^C F C A A 2 ■ 0.5314. (37) 

1X Z JO IX 

This correlator describes cross-talk between the two hemispheres and it is positive. As we will 
see below, this has direct consequences for the differential event shape distributions, confirmed by 
phenomenological analysis [0. 

So far, we have kept the couplings fixed in our investigation of {{£r£l)). The inclusion of higher 
order logarithmic corrections, of course, requires a resummation, and we certainly do not regard 
the numerical value of the correlation in ( |3TD as a prediction. The positive sign of ((£r£l)), 
however, is stable under any resummation that incorporates long-distance effects through the 
running of the couplings with transverse momenta, and even as functions of the relative rapidity 



rj. This follows from the form of integrand of Eq. (33), which is positive-definite for any fixed 
values of the transverse momenta and 7]. Similarly, the explicit Xrdependence in (f28|) is a direct 
consequence of boost invariance, and holds both perturbatively and nonperturbatively. 

5. From the correlators to the shape function. We are now ready to apply the expressions 
that we have obtained for the energy flow correlators to reconstruct the shape function f(en, £l)- 
We start with the single-gluon approximation, in which the factorization property (p6|) holds, and 
for which the correlations of energy flow are given by Eqs. (p4[) and (|25|). It follows from (p6[) that 

/OO rl\ 
7TZ lX£R (expHA^)) , (38) 
-00 2,71 

and similarly for /(^l)- To calculate the expectation value entering this expression, we use the 
cumulant expansion, 

00 1 

<e X >=exp£-«A^», (39) 
j= i3- 
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where the lowest cumulants are, ((1)) = (1), ((12)) = (12) - (1)(2), ((123)) = (123) - (1)(23) - 
(12)(3) — (13) (2) + 2(1)(2)(3), and analogously for higher cumulants (see, for example fL6|). In 
introducing the cumulant expansion, we note the close resemblance to the classic analysis of 
particle multiplicity, |jTT| . 

The shape function /(er) is found by substituting X = —iX£ R in (|39|). In the one-gluon 
approximation, ((£ N )) = (£ N )- Then, using the correlators (|24"D , we easily resum the series in 
(j38|) and obtain 

f(e R ) = jT ^ exp Ue r + jf dk 2 p PT (k 2 ) jT d V (e-^'" - l) j , (40) 

and analogously for /(el)- It is straightforward to verify that /(er^l) is a positive definite distri- 
bution. The function we have just calculated describes a resummed single-gluon contribution to 
the shape function. Its lowest-order expansion in powers of a s reproduces a familiar 1/er energy 
distribution, 

f(E R ) = —9 {A - er) f dk 2 p PT (fc 2 ) + O (a 2 ) . (41) 

At small er, higher-order corrections become important, and the asymptotic behavior of the 
resummed expression is drastically changed. To see this, we notice that the er — > behavior of 
(fiO|) is determined by the asymptotic behavior of the integrand as A —>■ oo. In this limit, using 
/ °° drj jexp (— iXy/k^e' 11 ^ — 1 j oc — In A, we find 

f{e R ) ER ~° eT~\ ap T = £ 2 dk 2 p PT (k 2 ) . (42) 

Here opt has the interpretation of the number of particles per unit rapidity interval. Perturbation 
theory thus suggests that / vanishes as a power of the energy, although it cannot reliably predict 
the absolute value of the exponent opt- In particular, apT suffers from infrared renormalon 
ambiguities related to the divergence of the coupling constant at small values of k 2 . One can 
estimate onp ~ 1 — 3, replacing p PT by a phenomenologically motivated pnp P"3"| . 

Now let us consider slightly larger values of Er, for which Eq. (fiOl) is dominated by finite values 
of A, that is, values such that X\/k? is a number of order unity for all k 2 < A 2 . In this region, we 
may expand the exponential, exp (— i\\^k^e^ v ^j as a power series in A. If we keep only the linear 
term in A, we get 

fW(e R )o<L6(e R -(£ R )) , (43) 

where the superscript refers to the linear approximation used in (flOD. This model for the shape 
function is equivalent to a "shift" of the perturbative distribution in Eq. (|T0l). 

Clearly, a simple shift cannot readily be combined with the \ow-er behavior of Eq. (f42"D. On 
the other hand, the shift results from keeping only the lowest order in A in our expansion of the 
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exponent in Eq. (fJOJ). The inclusion of the next order, A 2 , term amounts to a smearing of the 
delta-function into a Gaussian form 

« , ^— exp {-^4vr} , e R = - (44) 

where ((£#)) = + C^a 2 ), with the first term given by fl2lj) . 

At order A 2 , we may also begin to take into account the effects of the correlations between 
hemispheres, Eq. (pT7|), in the shape function (^). Using the integral representation for the delta- 
function entering (||), and applying the expansion ([39]) we get, for e R;L ~ A, 

/(e*,e L ) w [°° ^ ^ e^+W) exp (_ z((Al £ l + \ R S R )) - \{{{X L £ L + \r£r) 2 ))} , 

(45) 

where we have neglected contributions of higher cumulants. In this region, we find the approximate 
behavior 

f(e R ,e L ) « -LJ^expj-^ (<<£ 2 ))4 + <<£ 2 )) £ -| - 2((£ r £ l ))e r e l )] . (46) 
Here the dispersion is defined as 

V={(£l)){(£l))-((E R £ L ))\ (47) 

and clearly ({£%)) = ((££))■ For small £ R> l, the correlator ((£ R £ L )) does not affect the asymptotic 
behavior of the shape function. Using (f4~2|) one finds f(e R ,ei,) ~ /(£r)/(£l) ~ (£R£i) aNP ~ as 
£.r,l "~ > 0. Together, these results imply, almost uniquely, the parameterization, Eq. (|i3|) for the 
double-hemisphere shape function f(e R , El), as a minimal form, neglecting higher cumulants. The 
parameter a of (|T3| ) has the interpretation of the particle multiplicity per unit rapidity radiated 
by the boost-invariant sources, while b measures the extent to which radiation in one hemisphere, 
produced perturbatively or otherwise, "spills over" into the adjoining hemisphere, increasing cor- 
relations between e r and Ej_,. 

The Gaussian fall-off of (|13|), of course, need not extend to arbitrarily large values of El,r- 
For example, the large-e^ behavior that follows from the one-gluon approximation Eq. (fiO|) is 
actually ~ exp{— (e/A) \b.(e/A)}, which decreases somewhat less rapidly than the Gaussian in 
Eq. fll3|), although faster than an exponential. From a phenomenological point of view, however, 
this difference should have a modest effect on convolutions like Eq. flTU| ) for physical cross sections. 

To estimate the effect produced on the shape function fl4~6|) by nonzero correlations between 
the right and left hemispheres, we examine two extreme cases, ((£ r £l)) = ((££)) an d ((£r£l)) = 
— corresponding to maximal and minimal value of the correlator, respectively. Since in both 

cases the dispersion ([IT]) vanishes, the shape function is reduced to S(er— El) and 5(e r +El— 2(£ r )), 
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correspondingly, indicating that the energy flows into the two hemispheres are strongly correlated 
with each other. For ((£r£l)) > this correlation is positive - as Er increases, the same is true 
for el- For ((£r£l)) < the situation is opposite - the energy flows into two hemispheres are 
anti-correlated with each other. 

It follows from our calculation, Eq. fl3~7|), that the correlations between two hemispheres are 
positive. This property is in agreement with phenomenological fits of a Gaussian-type ansatz to 
experimental data || on event shapes. It is also in accord with the dispersive approach |18 . 



6. Conclusions. The event shape functions acquire an interpretation of energy flow distri- 
butions through hemispheres at infinities separated by the plane orthogonal to the jet axis. We 
have used perturbative QCD as a tool in our analysis of these functions, supplemented by con- 
ventional assumptions on the correspondence between perturbative and nonperturbative behavior 
||. While the resulting rapidity dependence is realistic, the transverse momentum dependence is 
governed by nonperturbative physics. The resumation of multi-energy correlators, in a one-gluon 
approximation, results in a power fall-off as Er^l — > 0. The sign of the correlation between the left 
and right hemispheres is more general than its perturbative approximation, and is a manifestation 
of the tendency of radiation into one hemisphere to spread into the other. As such, it is sensitive to 
the hadronization dynamics of soft radiation at wide angles from the jets. The functional form of 
the shape functions resulting from our approach is in the agreement with recent phenomenological 
analyses of differential event shape distributions || [7|. We expect to extend our results to all 
orders in the underlying perturbative expansion. 
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